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^ \ Abstract 

Following earlier one-step calculations, we explore the contributions of multi- 
step effects for the breakup of low energy *B on ^^Ni and ^'^^Pb within a 
coupled discretised continuum channels (CDCC) formalism. The Coulomb 
ps| , multi-step differential cross section is significantly reduced for all angles, the 

^^ \ largest effect being the destructive interference of nuclear couplings. The 

^) ' nuclear peak, at around 80° in the one-step calculations for ^^Ni, virtually 

On ■ disappears. 

0\ ■ 

_^ 24.10.Eq, 25.60.-t, 25.60.Gc, 27.20.+n 

-(— > ■ 

> 

X 



Typeset using REVTeX 



*Email address: filomena@wotan.ist.utl.pt 

1 



I. INTRODUCTION 

Coulomb dissociation has been proposed on many occasions |I|J^ as a means of determin- 
ing the interaction between fragments at low relative energies. This method is applicable 
even when the fragments are themselves radioactive and not easily produced as targets for 
direct scattering experiments. It was believed that at sufficiently forward angles, and/or 
at sufficiently low energies, the impact parameters for the breakup trajectory would be 
large enough for Coulomb mechanisms to dominate, and for first-order theories of Coulomb 
breakup to be adequate. The Coulomb dissociation method deduces the radiative capture 
cross section by measuring the reverse reaction, the dissociation of a projectile (the fused 
system) by the Coulomb field of a target. 

The Coulomb dissociation method has been used to examine the breakup of ®B at both 
high p,Q and low [^] energies. Analyses have started with semiclassical theory 0], and have 
progressed to include the E1/E2/M1 contributions with correct experimental efficiencies 
1^, three-body kinematics in the final state [^, and most recently one-step nuclear and 
Coulomb contributions ||9HTI|. These last results, and those of ref. [|12|, showed that the 



nuclear and Coulomb form factors extend to considerably larger distances than the sum of 
the radii of the participating nuclei, because of the extended tail of the wave function of 
the last proton in ^B, and that there is strong Coulomb-nuclear destructive interference at 
intermediate radii. We are prompted, by the size of these effects, to examine the importance 
of multi-step contributions for both Coulomb and nuclear processes, taking into account the 
final state interactions which were previously omitted. These final-state interactions will 
couple together the different continuum states, and also describe the depletion of the elastic 



channel due to breakup. The depletion effect has been considered in some calculations [[12 
where only couplings between the bound state and continuum states were included. The 
contribution of 'higher-order breakup' has yet to be properly clarified, especially the role of 
continuum-to-continuum couplings. 

Earlier treatments of multi-step effects for breakup have used either adiabatic [|13| or 



semiclassical approximations |14,15|, solved the scattering problem with time-dependent 



methods [0, tried CDCC solutions |T^JT^, or used Bremstrahlung integrals [|I^,^. Esti- 



mates of second-order Coulomb and nuclear effects have also been calculated [pT |. 



The earliest adiabatic approximation used the three-body wave functions of Amakawa 



et al. [|22| within a prior-form breakup matrix element [^ for both nuclear and Coulomb 
mechanisms in the ^Li breakup into a + t, when incident on 208p|^ ^^ jq MeV. However, 
the best fit to the experimental cross sections was found when both Coulomb distortion and 
Coulomb breakup were omitted from the calculation. The breakup of ^Li into a + d, when 
incident on 208p|^ ^^ ]^gg MeV, was later analysed [|T^ within the CDCC framework []TB 



and strong nuclear and Coulomb interference effects were found, even at forward angles. We 
follow a similar CDCC approach, but extend it to include, for the first time, dipole as well 
as quadrupole Coulomb mechanisms. 

Recent investigations of the breakup of halo nuclei have prompted a revival of semiclas- 
sical treatments of breakup p!^ , p!5| , p3| where the continuum is discretised into an orthogonal 
set of basis functions. First-order and higher-order couplings can then be included when 
integrating along a semiclassical trajectory (Rutherford orbit at low energies, or straight 
lines at high energies). Simplified ground-state wave functions are often used, and collective 



rather than semi-microscopic form factors calculated. We will see below that both of these 
approximations have to be reviewed when we consider the breakup of ^B incident on ^^Ni 
at the sub- Coulomb energy of 26 MeV. 

Another method, successfully used [p!6| , ^ to treat higher-order processes, is to follow the 



breakup reaction as a time-dependent process. Esbensen et al. |T6| follow the time evolution 
of a ^Be+p bound state by means of a TDHF propagator along a straight-line trajectory, 
and can calculate both Coulomb and nuclear contributions to breakup in a unitary manner. 



The results [0 show that the Coulomb higher-order dynamical processes cause a destructive 



E1/E2 interference, and a reduction of the dissociation probability. 



Most recently, the adiabatic three-body wave functions have been used again 1 19,20 



for Coulomb breakup, since in this case analytic solutions have been discovered both for 



the three-body wave functions |^ and for the post-form T-matrix integral of the breakup 
matrix element |]19| in terms of a Bremstrahlung integral. Unfortunately, the method is not 
immediately applicable for ^B breakup, since these analytic solutions only hold for neutral 
valence particles, and the adiabatic approximations are for high energy rather than for 
sub- Coulomb reactions. 

We have been progressively improving our understanding of low energy breakup reactions 
P|^TT|, focusing in particular on the breakup of ^B on ^^Ni measured by the Notre Dame 
group at 26 MeV. The work presented here is an exploratory continuation of our previous 
investigations, and now, for the first time, a full multi-step quantum mechanical description 
is attempted, including all continuum couplings. 

II. THEORY 

A. Coupled Discretised Continuum Channels (CDCC) 

When a projectile is described as a single particle outside a core, its state can be disturbed 
by the interaction with the target nucleus, as the tidal forces of the target act differentially 
on the particle and the core. If one separates the projectile-target interaction into V^t(R-c)) 
the interaction of the target nucleus with the core, and Vpt{r'), the interaction of the target 
nucleus with the particle, then there is a mechanism for coupling ground and inelastic 
(continuum) states together. Nuclear and Coulomb components of Vd and Vpt should be 
included on the same footing. 

In order to describe the breakup of a projectile such as ^B, we could consider the inelas- 
tic excitations in the p+'^Be system from the ground state (pgs{r) to excited states in the 
continuum Uesj,k{f'), for some momentum k and partial wave i. The use of such single energy 
eigenstates, however, would result in calculations of the inelastic form factors which will not 
converge, as the continuum wave functions do not decay to zero as r — * oo sufficiently fast to 
be square integrable. One way PBlETf of dealing with this divergence is to take continuum 



states, not at a single energy, but averaged over a narrow range of energies, such that these 
'bin' states are square integrable. We label these bin states by their wave-number limits 
[ki, ^2] and their angular momentum quantum numbers {is)j. We use them in the coupled 
discretised continuum channels (CDCC) method [^,^. 

Let R be the coordinate from the target to the projectile and r the internal coordinate 
of the projectile. The position coordinates of the projectile fragments with respect to the 



target are 



and 



R+^^ 



IXr — R. 



Ar. 



Ar. 



where the projectile has mass Ap. 

The wave function for the three-body system of proton + ^Be core + target is expanded 
as 



^-(R,r) 



1 






with 
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:r) = -J2{imsa\jf,)Yr{r)x:u, 
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V) 
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where Xs is the proton's intrinsic state. The set {L, £, s,j, [ki, ^2]} will be abbreviated as a. 
The radial wave functions fL,ij[ki,k2],j{R) satisfy the set of coupled equations 



h' 



<f L{L + r 



2/i V^i?2 



i?2 



+ ei[h,k2])-E 



f^j{R) + J2t^-^ ViAR)U,j{R) = 0, (3) 



where e{\ki,k2\) is the average energy of continuum bin [A;i,/c2], (or e < for the ground 
state). V^.^i{R) describes the coupling between the different relative motion states: 



K.;a' = < </'a(r)|l4t(Re) + Vpt{v')\ct>o.,{v) >, 



(4) 



where V^j(Rc) and Vptij') are again the total (nuclear and Coulomb) interactions between 
c — t and p — t systems respectively. In Eq. (^ radial integrations are done over r from zero 
to Run^ a parameter to be chosen. 

The coupled equations of Eq.(|) may be solved exactly |^ if they are not too numerous. 
Otherwise, iterative expansions are used starting with f^~^\R) = 0, and continuing as 



n' 
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for n=0, 1, ... The function /'•°^(i?) is thus the elastic channel, and the asymptotic S- 
matrix S^"''^ of the wave functions /'■"^(i?) gives the cross section for n*''-order DWBA. The 



n = 1 first-order DWBA solutions are presented in the previous paper |]Tl|. The multi-step 
DWBA results for large n will converge to the coupled-channels solution if the off-diagonal 
couplings are small. If they are large, the DWBA series diverges. Then the infinite series 
may be summed by the method of Fade approximants described below. 
The bin wave functions are defined as 



Uisj,[ki,k2]i'^) = \j ^ / w{k)e''^''uesj,k{r)dk , (6) 

with 6k the scattering phase shift for Uisj^ki^), the single-energy scattering wave function 
in the chosen potential Vp^{r) which may be ^-dependent. The normalisation constant is 
N = Jf^^ \w{k)\'^dk for the assumed weight function w{k), here taken to be unity. These bin 
states are normalised {u\u) = 1 once a sufficiently large maximum radius Rbm for r is taken. 
They are orthogonal to any bound states, and are orthogonal to other bin states if their 
energy ranges do not overlap. The phase factor e~*'''= ensures that they are all real valued 
for real potentials Vp^{r). 

The rms radius of a bin wave function increases as the bin width ^2 — ki decreases, 
approximately as 1/(^2 — ^i), so large radial ranges are needed to include narrow bin states. 
If the maximum radius Rbm is not sufficiently large, then the bin wave functions U[ki,k2] will 
not accurately be normalised to unity by the factors given in equation (^). It is important 
however, to realise that the missing normalisation comes at large distances; the bin wave 
functions must not be artificially renormalised to unity, otherwise, for example, the correct 
Coulomb B{EK) distributions will not be obtained. 

The couplings V^.^,{R) in Eq.(||) arise, as discussed above, from the interaction potentials 
of the projectile fragments with the target. Assuming that the potentials Vd and Vpt are 
central, the Legendre multipole potentials can be formed as 

Vk{R, r) = ^ J^' [V.ti'Rc) + Vptir')] PK{x)dx , (7) 

where K is the multipole and x = f ■ R is the cosine of the angle between r and R. Since 
s is the (fixed) spin of the proton, a spectator, the coupling form factor between states 
Ui'ikiM'i'^) ^^^ ^i[ki,k2]i'^) is then 

VL'iR)=E i-iy^''-'-'jj'ii'LL'i2K + l)W{jj'ii'; Ks)W{jfLL'- KJ) 



K 



K i i'\ f K L L'\ /•«"- 



^lo Oojio I Jo '^isjlk^Mli^) ^K{R,r) Ue^sj'lkiMvir) dr ■ (§) 

From the S'-matrices S**-"-* we calculate the double differential cross sections d^a /dVtde^ 
where VL is the scattering angle of the centre of mass of the ^B* fragments p and c, and e is 
the excitation energy within ^B*. Usually we will plot the integrated angular distribution 
da/dfl obtained after summation over all the bin energies e([/ci, ^2])- It will be possible 
Il27| , |28[| to obtain from the CDCC results the full multiple differential cross sections such as 
d^a / dQpdflcdEp, and then, for example, any post-acceleration effects may be determined 
from the variation of the cross section with respect to the laboratory proton energy Ep. An 



approximate treatment is available [p,ll| to calculate these fragment distributions using an 



isotropic assumption for the pc relativ motion, but this does not yet include the interference 



between the different final ^-states that is necessary |T5[ to give non-zero post-acceleration 
effects. 



B. Fade acceleration 

A given sequence S^^\ S'^^\ ■ ■ ■ of S-matrix elements that result from iterating the coupled 
equations can be regarded as the successive partial sums of the polynomial 

/(A) = 5^ + (5(1) - 5W)A + (5(2) - 5(i))A2 + ■ ■ • (9) 

evaluated at A=l. This polynomial will clearly converge for A sufficiently small, but will 
necessarily diverge if the analytic continuation of the /(A) function has any pole or singular- 
ities inside the circle |A| < 1 in the complex A-plane. The problem that Pade approximants 
solve is that of finding a computable approximation to the analytic continuation of the /(A) 
function to A=l. This is accomplished by finding a rational approximation 

"^^^'^^^ ^^^ = l + g,A + g,A^ + ... + gMA-^ ^'°^ 

which agrees with the /(A) function in the region where the latter does converge, as tested 
by matching the coefficients in the polynomial expansion of P[n,m]{^) up to and including 
the coefficient of A" for n = A^ + M. 

There are many different ways |^| of evaluating the coefficients pi, qj, but for the present 



problem we can use Wynn's e-algorithm |31|j3^ , which is a method of evaluating the upper 
right half of the Pade table at A=l directly in terms of the original sequence S^^\ S^^\ ■ ■ ■. 
Experience has shown that for typical sequences the most accurate Pade approximants are 
those near the diagonal of the Pade table. We use S^") = P[n,m]{1) for N = [{n + l)/2] and 
M = [n/2] in calculating the Pade-resummed cross sections. 

When accelerating a vector oi S-matrix elements S'-"), with a component for each coupled 
channel a, then it is important to accelerate the vector as a whole. Wynn [^| pointed out 
that this can be done using the Samuelson inverse 

x-i = (x-x*)-ix* (11) 

where x* is the complex conjugate of x. 

III. RESULTS 

The multi-step DWBA calculations presented here were calculated using Fresco pU 



We use a continuum breakup subspace sufficient to reproduce what we believe to be the 
principal channels. For the distorted waves of the projectile-target wave function, radii 
up to Rcoup = 300 fm and partial waves up to Lmax = 600 were included to ensure full 
convergence of the individual angular distributions. These limits give, by semiclassical 
considerations, cutoffs for Coulomb excitations below 2.0° from our L^ax, and below 1.7° 
from our Rcoup value. We have examined the convergence with respect to Rbin- For the case 
we are interested in, the differential cross section remains unaltered as long as Rbin > 50 fm. 
Thus, each energy bin (Eq.| and Eq.|) is calculated using Rbm = 50 fm. 

It is essential for physical completeness that our calculations include monopole, dipole 
and quadrupole contributions for both nuclear and Coulomb mechanisms. However, they do 
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not include Ml transitions. At the extreme nonrelativistic velocities of interest here, these 
are predicted (see for instance p) to be insignificant. 

In Fig. (|I]) we show the energy distribution of the cross section obtained within a 1-step 
calculation, using the ^B model from Esbensen ^^. We keep the same p+^Be potential 



(that defined for the ground state) for all partial waves of relative motion. For scattering 
from ^^Ni we use the same optical potentials [p4| , ^ as in f^. The cross section is plotted 
as a histogram to illustrate the continuum discretisation that we have used to define the 
energy bins included in all calculations. These results show that transitions from the ^Be-p 
P3/2 ground state to s,p,d, and f-wave continuum states up to 3 MeV should be taken into 
account, even though one can expect f- waves to offer only a small correction to the overall 
result. A finer discretisation of p, d and f waves would be desirable, but the present (70 bin) 
calculation is already at the limit of our computational capacity. 

We first calculate multi-step effects by iterating the coupled equations beyond the first- 
order DWBA. We find that even the second-order DWBA diverges rapidly, especially for 
low partial waves (small impact parameters), and does not give sensible results beyond 30° 
(dashed line with circles in Fig. |^). In order, therefore, to present some indications on what 
may be deduced from the successive Born terms, we will present the results when resumming 
the expansion using the method of Fade acceleration (see section [11 Bp for successive numbers 



of steps n. We will use for cross sections the Fade approximants the S^") rather than the 
original S^"'^ matrices. Given the nature of our expansion, it is not possible to directly 
compare 2nd and 3rd order effects with those obtained using the pure conventional DWBA 
expansion. 

Fig.(0) shows the 1-step, 2-step, 3-step, 6-step and 20-step breakup results using Fade 
acceleration. The rate of convergence of this resummed expansion is encouraging, contrary 
to that of the original Born series which diverges strongly immediately at 2nd order. In 
addition, we show in figs.(|^) and (|^) the different rate of convergence for the ^B breakup 
into s and p continuum states considering the Coulomb and nuclear interactions separately. 
In all the cases we have studied, the Fade convergence is non monotonic. If one includes 1- 
step and 2-step processes only, the differential cross section is underestimated. Introducing 
3-step corrections overestimates the cross section. From our results we conclude that for 
the breakup of low energy ^B on ^^Ni, contributions up to at least 9th order in the Fade 
expansion should be included. 

Still in Fig.(^ we present the results obtained for the full coupled channel calculation 
taking into account s-wave and p-wave bins. We find that processes beyond 20-step do not 
contribute to the cross section. For this reduced bin subspace (42 bins) it is possible to 
perform the full coupled channel calculation (light solid curve in Fig.^. It is reassuring to 
find that our results using the multi-step expansion with Fade acceleration converge to the 
correct full CDCC results. In the larger bin subspaces it is extremely hard to perform the full 
CDCC calculation and thus we will rely on the multi-step expansion with Fade acceleration. 

Including multi-step effects, the Coulomb differential cross section is hardly modified up 
to 6 c^ 10°. The peak at ^ ^ 20° is shifted to slightly smaller angles with higher order 
processes and its magnitude is reduced by ~ 10% (see Fig.^. 

The most striking result of our work is clearly the destructive interference caused by the 
nuclear multi-step processes. The nuclear peak is shifted to lower angles (from 6 ~ 80° for 
1-step calculations to ^ ~ 40° for the CDCC calculations) and suffers a reduction to ^ of 



its peak value. We do not expect measurements of the breakup differential cross section at 
larger angles to provide a good handle for the optical potentials as one could deduce from 
the 1-step results presented in our earlier work |p. The previously observed strong nuclear 
peak is practically washed away by multi-step effects. 

We point out that, as in the 1-step calculations of 0, the total differential cross section 
does not correspond to the sum of Coulomb and nuclear contributions calculated separately. 
This can be seen in Fig.(^) where the sum of the Coulomb and nuclear cross sections for the 
CDCC calculation is compared with the CDCC cross section when Coulomb and nuclear 
are treated in the same footing. As in P,p!T| there is a wide range of angles where the 
Coulomb-nuclear interference effects cannot be neglected. 

So far we have included all possible couplings within the subspace considered. However 
it is useful to identify the relative importance of the continuum-continuum couplings as 
compared to the couplings to and from the ground state. In Fig.(||) we show the results 
of calculating the full multi-step breakup into s and p bin-states, including continuum- 
continuum couplings (dark lines) and excluding them (light lines). As can be seen, the 
continuum-continuum couplings are responsible for the significant cross section reduction, 
which are not merely due to depletion of flux from the elastic channel. The reduction is still 
obtained in calculations (not shown) with E2 couplings acting to only first order, but not in 
those with El couplings only to first order. This indicates that the reduction is caused by 
multistep El processes interfering with low-order E2 transitions, a process similar to that 



seen in ref. fl6 



In order to elucidate the different ^B partial wave contributions, we show in Fig.(^ 
the differential cross section obtained for the full multi-step breakup including: s (dotted), 
s-|-p (short-dashed), s+p+d (long-dashed) and s+p+d+f (dot-dashed) bin states. A good 
description of the physics can be obtained without d and f-waves, although if one wishes to 
extract quantitative results these should be included together with a finer energy-bin grid. 

One of the main motivations of ^B breakup experiments is astrophysical, to determine 
the 5*17 at low relative energies. It is thus important to disentangle the dependence on the 
'^B structure model. This was the main concern of our earlier work |jlO[. In Fig.(|^) we 
present a comparison of differential cross sections for two ^B models: that of Kim ||36|, and 



our initial model, from Esbensen [|T6l . The latter has a smaller radius in order to reproduce 
Sij ~ 17 eV b. The difference between the multi-step results using the two structure models 
is generally similar to the difference in the 1-step results: there is an overall normalisation 
due to the size of the projectile but no significant shape change. For this particular case, 
the Kim model produces a 30% increase in the total differential cross section. 

The Notre Dame group is considering repeating their experiment [Q with a heavier target 
at the same beam energy, hoping then the E2 contribution will be easier to extract. We have 
therefore performed one and multistep calculations for a ^osp^ target, for all combinations 
of nuclear and Coulomb mechanisms (see Fig.^). We have taken the ^ospb _ proton optical 
interaction from |^. For the optical potential between ^Be and ^°^Pb we have used a 



heavy- ion global parameterisation [g^. We have checked that the differential cross section 
is not sensitive to variations on the core-target optical parameters, and thus we expect these 
results to give a good indication of the physical effects. In order to have quantitative results, 
measurements of the elastic scattering of ^Be (or a nucleus in the same mass/charge region) 
on 208p]2 g^^ these low energies would be necessary. 



Our results show that the nuclear contribution is zero up to 50° and becomes important 
only at backward angles. Given these results, it should be possible to extract information 
on the magnitude of the electromagnetic components, as long as the detectors are placed at 
smaller angles. One should keep in mind that the multi-step processes reduce the Coulomb 
peak and alter its shape. 

Similarly to what was found for ^^Ni and 1-step DWBA, for a ^"^Pb target there are 
interference effects that do not allow a simple subtraction of the E2 component, as one 
would wish to obtain the 5*17. In Fig.(|l^) we show the results for the CDCC calculations for 
the different electric components together with the full calculation. One can clearly see a 
destructive interference between El and E2 components. It may be possible to disentangle 
these components, but the result will inevitably be model dependent. 

IV. CONCLUSIONS 

Multistep calculations of low energy ^B breakup on ^^Ni and ^°^Pb, including all relevant 
couplings, have been performed for the first time. We have found it necessary to treat nuclear 
and Coulomb potentials on the same footing, since there is considerable interference between 
these mechanisms. Our calculations use the CDCC method, with Fade approximants to 
resum the Born series for the S-matrix. Exact coupled-channel calculations, possible in a 
reduced subspace, verify that this Fade resummation converges to the full coupled channel 
result. We compare our results with the one-step prior- form DWBA cross sections reported 
previously ||^|ll[ . 



The multi-step effects are very strong, producing significant reductions of the cross sec- 
tion compared with those from first order theory. For Coulomb breakup we see pronounced 
interference effects for all angles at and beyond the peak position, while multi-step effects 
are much stronger for the nuclear part, so that, for the ^^Ni target, the nuclear peak resulting 
from the 1-step calculations virtually disappears. The prior-form DWBA thus overestimates 
the nuclear breakup probabilities at our sub-Coulomb incident energy. The multi-step re- 
duction is principally due to the continuum-continuum couplings, not just to depletion of 
the elastic channel, and this indicates that the projectile undergoes considerable dynamical 
distortion and recombination during the reaction. The dominant qualitative changes caused 
by the multi-step effects are seen when including only s and p wave continuum bins, but d 
and f waves must be included for quantitative results. 

The results for two different ^B structure models show that the multi-step effects depend 
on the size of the projectile, and are not sensitive to other details of the ®B g.s. wave 
function. With a ^osp]^ target, the nuclear contribution is only significant for backward 
angles. Multi-step effects reduce the cross section and change the shape slightly. According 
to our results, for the extraction of the E2 component this experiment seems more promising 
than that with the lighter target, although care should be taken to account for the strong 
destructive Coulomb- Coulomb interference. 
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FIG. 1. Energy distribution of s, p, d and f partial wave bins. 
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FIG. 2. The differential cross section obtained for multi-step breakup of '^B into s and p-wave 
bins, including both Coulomb and nuclear effects: the full CDCC calculation, the 1-step and 2-step 
DWBA, and higher order calculations using Fade acceleration. 
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FIG. 3. The differential cross section obtained for the multi-step Coulomb breakup of ^B into 
s and p-wave bins using Fade acceleration. 
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FIG. 4. The differential cross section obtained for the multi-step nuclear breakup of ®B into s 
and p-wave bins, using Fade acceleration. 
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FIG. 5. The CDCC differential cross section obtained for the breakup of ^B into s-wave and 
p-wave bins: comparison of the Coulomb and nuclear summed cross section with the calculation 
that includes nuclear-Coulomb interference. 
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FIG. 7. Cumulative ^B partial wave contributions to the CDCC differential cross section for 
Coulomb plus nuclear breakup to the 1=0,1,2 and 3 channels. 
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FIG. 8. Sensitivity to the ^B structure model: the CDCC differential cross section including 
all but f-wave bins. 
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FIG. 10. Comparing the CDCC effects of only EO, only El, and only E2 with the full Coulomb 
differential cross section for ^B Coulomb breakup on ■^''^Pb, including all but f-wave bins. 
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